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Introduction: 

The  primary  goal  of  this  project  is  to  develop  a  zebrafish  model  of  the  autosomal  dominant 
genetic  disorder  type  I  neurofibromatosis  (NF1).  This  disorder  is  very  common,  affecting 
approximately  1  in  3000  live  births,  and  results  from  mutations  in  the  NF1  gene.  A  zebrafish  model  of 
this  disease  will  be  particularly  useful  in  furthering  our  understanding  of  the  pathophysiology  of  this 
disorder  and  will  allow  the  application  of  high  throughput  strategies  to  screen  libraries  of  chemical 
compounds  to  identify  small  molecules  that  may  prove  efficacious  in  modulating  NF1 -associated 
phenotypes.  The  zebrafish  model  will  also  allow  for  rapid  and  cost-effective  structure-function  analysis 
of  neurofibromin,  the  large  protein  product  of  the  NF1  gene.  Zebrafish  develop  rapidly  and  are 
transparent  throughout  early  development  allowing  for  the  easy  visualization  of  many  developing 
tissues,  including  the  cardiovascular  and  nervous  systems.  In  addition,  the  availability  of  transgenic 
reporter  lines  that  express  fluorescent  proteins  under  the  control  of  tissue-specific  promoters  allows  for 
real-time  visualization  of  complex  developmental  processes.  This  model  system  is  also  tenable  to 
forward  genetic  screens,  which  will  allow  for  the  identification  of  enhancers  and  suppressors  of  NF1 
phenotypes  in  an  unbiased  fashion.  Prior  to  the  efforts  funded  by  this  award,  the  zebrafish  orthologues 
of  the  human  NF1  gene  had  not  been  described  or  characterized.  This  project  involves  the  collaboration 
of  two  established  groups  to  coordinate  efforts  in  developing  and  exploiting  a  zebrafish  model  of  this 
important  genetic  disease. 


Body  of  work: 

The  coordinated  efforts  of  the  Look  and  Epstein  laboratories  have  led  to  significant  inroads  into 
the  primary  objectives  of  this  award.  We  outline  below  our  collective  progress  since  inception  of  this 
work  reiterating  aspects  of  our  investigations  described  in  our  last  annual  report  that  have  been 
significantly  expanded  upon  along  with  new  information  not  previously  presented. 

Identification  of  nfla  and  nflb 

We  utilized  a  bioinformatics  approach  to  identify  the  zebrafish  orthologues  of  human  NFL 
Analysis  of  the  eighth  assembly  (Zv8)  of  the  zebrafish  genome  revealed  two  genes  highly  similar  to  NF1 
at  the  amino  acid  level  (90.4%  and  90.7%  respectively),  which  we  named  nfla  and  nflb.  These  genes 
are  highly  related  to  one  another  (87.4%  identical  and  93.7%  similar)  with  nfla  and  nflb  sharing  similar 
genomic  structures  and  each  containing  57  exons  (Fig  1A-B).  nfla  is  located  on  chromosome  15  (Fig 
1A)  and  predicts  a  311  kD  protein  composed  of  2755  amino  acids  while  nflb  is  located  on  chromosome 
10  (Fig  IB)  and  predicts  a  310  kD  protein  composed  of  2747  amino  acids. 

Comparison  to  Drosophila,  murine,  and  human  neurofibromin  protein  sequences  reveals 
significant  conservation  in  the  GAP  and  IRA  homology  domains  and  also  in  extensive  areas  flanking 
these  regions,  suggesting  additional  functional  motifs  that  have  been  conserved  across  evolution  (Fig  2). 
A  phylogenetic  tree  (Fig  1C)  demonstrates  a  tight  clustering  of  the  zebrafish  neurofibromin  orthologues 
with  other  mammalian  neurofibromins  and  a  divergence  from  the  Drosophila  neurofibromin  orthologue. 
Human/zebrafish  synteny  maps  and  bioinformatics  analyses  suggest  that  nfla  and  nflb  likely  arose  via 
gene  duplication  (Fig  ID).  Upstream  of  the  human  NF1  gene  on  chromosome  17  are  genes  encoding 
WD  repeat  and  SOCS  box-containing  1  (WSB1),  Kinase  suppressor  of  Ras  1  (KSR1),  and  Galectin-9 
(LGALS9)  while  A  kinase  anchor  protein  1  (AKAP1)  and  RNA-binding  protein  Musashi  homolog  2 
(MSI2)  both  lie  downstream  of  NF1.  Similar  genes  flank  nfla,  while  nflb  is  flanked  only  by  orthologues 
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of  KSR1  and  MSI2.  The  identification  of  duplicated  genes  is  common  in  zebrafish  and  reflects  the  well- 
described  chromosomal  doubling  event  occurring  early  in  teleost  evolution. 

nfla  and  nflb  are  expressed  maternally  and  in  the  developing  zebrafish  cardiovascular  system. 

The  GenBank  EST  database  identified  expressed  sequence  tags  for  both  zebrafish  nfl  genes  in 
many  tissues  including  the  heart  (Fig  4A)  suggesting  that  neither  gene  is  likely  to  be  a  pseudogene  and 
that  they  are  expressed  in  overlapping  tissues.  We  examined  the  expression  of  both  genes  by  whole 
mount  in  situ  hybridization  between  the  4-cell  stage  and  4  days  post  fertilization  (dpf)  and  found  that 
both  genes  are  expressed  ubiquitously  during  early  development  with  later  restriction  to  regions  of  the 
head  and  anterior  central  nervous  system  (Fig  3A1-7,  Bl-7,  Fig  S2C1-C7,  D1-D7).  Notably,  at  48  hours 
post  fertilization  (hpf)  and  3  dpf  both  genes  are  expressed  in  the  heart  (Fig  3  A4,  B4,  A6,  B6,  Fig  4C5-6, 
D5-6)  and  in  the  dorsal  vessel  (Fig  3A3,  B3,  A5,  B5,  Fig  4C7,  D7).  Reverse  transcription  and 
polymerase  chain  reaction  (RT-PCR)  using  RNA  from  wild  type  24,  72,  and  84  hpf  whole  embryos  or  3 
dpf  Tg(kdrl:GRCFP)znl  GFP-positive  sorted  cells  confirmed  expression,  particularly  in  the  vascular 
endothelium  (Figure  4E-G),  while  RNA  from  1-cell  embryos  indicate  that  both  genes  are  expressed 
maternally  (Fig  4B).  Queries  of  an  expression  database  generated  from  sorted  endothelial  cells  from 
Tg(flil:egfp)yl  zebrafish  identifies  nfla  and  nflb  in  both  GFP+  and  GFP-  cell  populations,  consistent 
with  the  expression  of  these  genes  in  vascular  endothelium. 

Morpholino  knockdown  of  nfla  and  nflb  results  in  cardiovascular  defects 

nfla  and  nflb  morphant  embryos  displayed  gross  abnormalities  of  cardiovascular  development 
appreciable  to  the  blinded  observer  by  48  hpf.  Frequently,  blood  was  seen  to  move  back  and  forth  from 
atrium  to  ventricle  in  morphants,  suggesting  a  malfunctioning  atrio-ventricular  valve.  At  the  resolution 
afforded  to  us  by  histological  analysis,  we  observed  no  readily  apparent  structural  defects  in  the  atrio¬ 
ventricular  valves  of  morphants  despite  the  observed  functional  deficits  (Fig  5).  In  addition,  we  observed 
pooling  of  blood  in  the  common  cardinal  vein  and  a  paucity  of  blood  flow  along  the  dorsal  aorta  and 
posterior  cardinal  vein.  Valvular  insufficiency  and  reduced  blood  flow  were  not  seen  in  control 
morphants  or  wild  type  embryos.  Overall  development  of  the  embryos  was  relatively  preserved  through 
the  first  3  days  despite  these  cardiac  defects.  Histological  analysis  revealed  a  thinned  ventricular 
myocardium  and  large  pericardial  effusions  in  MO-treated  embryos  (Fig  6K-F,  Fig  7A-B,  F-G). 
Immunohistochemical  analysis  of  3.5  dpf  nfla,  nflb,  and  nfla/nflb  morphant  zebrafish  also  reveals 
increases  in  phospho-Erk  staining  (Fig  6M-N,  Fig  7C-E,  H-K).  Gross  morphological  analysis  showed  an 
increased  incidence  of  pericardial  effusions  beginning  at  48  hpf,  reflecting  cardiac  dysfunction,  in  nfla 
and  nflb  morhpants  when  compared  with  controls  (Fig  6A-G,  Fig  7L-M).  Non-specific  toxicity  arising 
from  MO  exposure  as  a  cause  of  the  observed  cardiovascular  defects  was  unlikely  as  unrelated  control 
or  scrambled  MOs  failed  to  produce  similar  levels  of  abnormalities,  defects  were  observed  even  at  low 
doses  of  specific  MOs,  and  similar  defects  were  observed  with  several  unrelated  but  specific  MOs 
directed  against  nfla  and  nflb.  In  addition,  injection  of  specific  MOs  in  p53  mutant  embryos  also 
produced  similar  cardiovascular  defects  (Fig  6H-J),  and  off-target  effects  due  to  MO  exposure  are 
known  to  be  partially  mediated  through  p53  activation.  Defects  in  cardiac  valve  morphogenesis  and  a 
thinning  of  the  ventricular  myocardium  are  also  seen  in  Nfl~'~  murine  embryos. 

We  performed  knockdown  experiments  using  zebrafish  embryos  in  which  endothelial  cells  are 
marked  by  expression  of  a  cytoplasmic  enhanced  green  fluorescent  protein  (GFP)  in  order  to  allow  for  a 
more  detailed  analysis  of  vascular  development.  Dramatic  abnormalities  of  vascular  patterning  in  the 
intersomitic  vessels  of  morphant  embryos  were  seen  at  48  and  72  hpf  (Fig  8).  In  nfla  MO-treated 
embryos,  the  leading  edge  of  the  sprouting  vessels  displayed  claw-like  projections  at  48  hpf  (Fig  8C) 
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and  failed  to  pattern  normally  such  that  the  dorsal  longitudinal  anastomotic  vessel  (DLAV)  did  not  form 
or  developed  in  a  rudimentary  fashion  (Fig  8F).  This  occurred  in  embryos  that  were  otherwise  normal  in 
overall  size  and  maturity.  These  defects  were  also  noted  in  nfla/nflb  compound  morphants,  and  were 
present  but  less  severe  in  nflb  morphants.  Vascular  patterning  defects  did  not  appear  to  correlate  directly 
with  cardiac  defects,  as  we  observed  embryos  with  vascular  abnormalities  that  did  not  display 
pericardial  effusion  or  valvular  insufficiency  as  assessed  by  a  “to-and-fro”  movement  of  blood  within 
the  heart.  Blood  flow  within  the  dorsal  aorta  and  posterior  cardinal  vein  appeared  intact  in  these 
embryos. 

At  24  hpf,  analysis  using  zebrafish  embryos  expressing  a  nuclear-localized  GFP  in  endothelial 
cells  indicated  that  morphants  displayed  a  complete  (Fig  9A4)  or  partial  absence  of  intersomitic  vessels 
emanating  from  the  dorsal  aorta  when  compared  with  stage-matched  controls  (Fig  9A3).  Overall 
morphology  of  morphant  and  control  embryos  appeared  equivalent  (Fig  9A1-2)  ruling  out  non-specific 
developmental  delay.  These  defects  were  apparent  following  MO-mediated  knockdown  of  nfla  or  nflb 
while  knockdown  of  both  together  had  an  additive  effect  (Fig  9B).  The  small  percentage  of  embryos 
with  defects  produced  by  the  5MP  MO  may  have  been  due  to  low-level  knockdown  of  nfla.  Again, 
similar  defects  were  observed  with  several  unrelated  but  specific  MOs  directed  against  nfla  and  nflb 
(Fig  10A).  Our  analysis  of  morphant  embryos  at  24  hpf  also  revealed  a  caudal  vessel  defect.  Morphant 
embryos  displayed  a  cystic  expansion  in  the  region  of  the  caudal  vein  and  exhibited  inappropriate 
anastomoses  between  the  caudal  vein  and  artery  (Fig  10D2-D4)  when  compared  with  controls  (Fig 
10D1).  Identity  of  the  expanded  tissue  as  vascular  was  confirmed  by  expression  of  GFP  (Fig  10D6- 
ZFN_MUT2)  and  the  observation  of  a  pooling  of  red  blood  cells  in  the  expanded  region.  This  defect  was 
present  following  knockdown  of  nfla,  nflb,  or  both  together  (Fig  10C-D). 

Additional  confirmation  of  the  role  of  nfla/nflb  in  vascular  development  derives  from  studies 
using  a  genetic  background  sensitized  to  vascular  insult.  Previous  studies  employed  MOs  directed 
against  flt4,  the  zebrafish  VEGF  receptor-3  orthologue,  to  investigate  genetic  interactions  during 
zebrafish  artery  development.  Additionally,/?^  morphant  zebrafish  embryos  display  variable  defects  in 
segmental  artery  formation  reminiscent  of  those  identified  in  our  nfla/nflb  morphants.  Endothelial-GFP 
expressing  zebrafish  embryos  were  injected  with  flt4  MO  alone  and  in  combination  with  a  MO  directed 
against  nfla,  nflb,  or  a  combination  of  both.  At  low  MO  doses,  85%  of  flt4/nfla,  24%  of  flt4/nflb,  and 
36%  of  flt4/nfla  +  nflb  compound  morphants,  displayed  abnormal  vascular  shunts  at  48  hpf  compared 
to  only  3-8%  of  individual /?r4,  nfla,  or  nflb  morphants  (Fig  9C,  Fig  10B).  This  defect  was  not  apparent 
in  control  morphants.  The  shunts  occur  between  the  dorsal  aorta  and  the  dorsal  longitudinal  anastomotic 
vessel  with  retrograde  flow  through  segmental  arteries  back  into  the  dorsal  aorta  or  through 
intersegmental  veins  into  the  posterior  cardinal  vein.  In  some  cases,  there  were  interruptions  of  the 
dorsal  aorta. 

nfl  loss  promotes  OPC  migration. 

We  have  shown  previously  that  oligodendrocyte  precursor  cell  (OPC)  numbers  are  increased 
upon  transient  nfl  knockdown  by  morpholino.  Here  we  report  our  analysis  of  the  behavior  of  OPCs 
following  nfl  knockdown  by  performing  live  time-lapse  imaging  of  these  cells  in  vivo.  The  oligl-EGFP 
transgenic  animals  were  monitored  for  12  hours,  beginning  at  60  hpf  when  OPCs  start  to  actively 
migrate  away  from  the  ventral  spinal  cord.  Compared  to,  control  p53  e7/e7  animals,  more  OPCs 
migrated  into  the  dorsal  spinal  cord  in  nfla+lb- knockdown  p53  e7/e7  embryos  (Fig  11  A,  B),  consistent 
with  our  observations  from  fixed  sections  (see  our  2009  progress  report).  Interestingly,  the  time-lapse 
study  also  revealed  that  nfla+lb- morphant  OPCs  traveled  longer  distances  relative  to  control  animals 
(Fig  1  IE).  When  representative  migratory  traces  of  individual  OPC  cells  were  analyzed  from  the  time- 
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lapse  movies,  nfla+lb- knockdown  OPCs  were  found  to  migrate  further  in  the  dorsal  and  rostro-caudal 
directions  than  control  OPCs  (Fig  1 1C  and  D).  Analysis  of  individual  OPC  cell  movements  over  the  12- 
hour  imaging  period  showed  a  35%  increase  in  the  total  distance  traveled  by  the  progenitors  in  nfl 
knockdown  animals  (32.05  pm  in  control  vs.  43.27  pm  in  nfl a+1  b-morphant,  p<0.05;  n=9  and  n=23, 
respectively;  Fig  11E).  In  all  cases  the  OPCs  display  intermittent  movement  consisting  of  repeated 
cycles  of  active  migration,  separated  by  stationary  pauses  that  often  distinguish  new  trajectories  from  the 
preceding  path.  After  nfl  knockdown,  the  OPCs  paused  for  shorter  periods  of  time  relative  to  controls 
(413.9  min  in  control  vs.  324.3  min  in  nfla+lb- morphant;  p<0.005;  Fig  1  IF),  however,  the  frequency  of 
the  pauses  did  not  differ  significantly  (2.28/hr  in  control  vs.  2.58/hr  in  nfl  a + 7/>morphant;  p>0.1). 
Furthermore,  the  velocities  of  OPCs  during  active  migration  were  determined  and  found  to  be 
unaffected  by  nfl  knockdown  (0.5203  pm/5min  in  control  vs.  0.5343  pm/5min  in  nfla+lb- morphant; 
p>0.7;  Fig  11G).  Taken  together,  these  findings  demonstrate  a  novel  in  vivo  phenotype  of  OPCs  in 
nfl  a+1  b-morphants  where  they  spend  more  time  actively  migrating  as  a  result  of  decreases  in  their 
period  of  inactivity. 

The  increase  of  OPC  number  upon  nfl  knockdown  depends  on  the  GAP  activity  of  NF1 

The  best-known  function  of  human  NF1  is  to  act  as  a  RAS-GTPase- activating  protein  (GAP)  to 
downregulate  RAS  signaling.  To  investigate  whether  the  increase  in  OPCs  seen  in  our  model  depends 
upon  the  GAP  activity  of  NF1,  we  performed  rescue  experiments,  injecting  mRNA  encoding  the  GRD 
of  human  NF1  into  nfla+lb- knockdown  embryos  at  the  one-cell  stage.  To  confirm  and  monitor  GRD 
expression,  we  fused  the  coding  sequence  of  the  GRD  domain  in-frame  with  the  gene  encoding 
mCherryRed  fluorescent  protein  at  its  N  terminus.  The  nfla+lb- knockdown  embryos  coinjected  with 
mCherryRed  RNA  as  a  control  showed  the  expected  increase  in  OPC  numbers  relative  to  control  MO- 
and  control  RNA-injected  embryos  (Fig  12B-D).  nfla+lb- knockdown  embryos  coinjected  with  the 
GRD  RNA  showed  a  relative  reduction  in  the  number  of  OPCs  in  the  dorsal  spinal  cord  from 
67.7/embryo  to  48.5/embryo  (n=10,  p<0.001;  Fig  12A,  B,  E).  Furthermore,  the  number  of  OPCs  present 
in  nfla+lb- knockdown  embryos  as  a  result  of  GRD  rescue  did  not  differ  significantly  from  the  numbers 
of  OPCs  present  in  control  MO-  and  control  RNA-injected  embryos  (48.5  vs.  47.5,  p>0.7,  and  48.5  vs. 
57.5,  P>0.1,  respectively).  These  results  indicate  that  the  increase  in  OPC  numbers  due  to  nfl 
deficiency  is  mediated  by  the  loss  of  GAP  activity. 

Generating  zebrafish  nfl  a  and  nflb  mutants  by  TILLING 

An  important  aspect  of  our  proposal  requires  the  generation  of  stable  zebrafish  lines  harboring 
mutations  in  nfl  a  and  nflb.  We  will  use  these  animals  to  verify  the  data  we  collected  from  our 
morpholino  studies  as  well  as  to  extend  our  analysis  to  time  points  beyond  those  accessible  using 
transient  knockdown  assays.  In  order  to  generate  stable  zebrafish  lines  carrying  mutations  in  nfl  a  and 
nflb  we  have  used  a  reverse  genetics  approach  called  TILLING  (Target  Induced  Local  Lesions  In 
Genomes),  in  collaboration  with  the  laboratory  of  Dr.  Solnica-Krezel,  and  have  found  several  nfl  a 
missense  mutant  alleles  and  one  nonsense  mutant  line  in  the  nfla  gene  (Fig  13A-B).  A  second 
mutagenesis  strategy  based  on  retroviral  insertion  identified  a  single  mutant  line  affecting  the  nfl b  locus 
in  the  intron  between  exon  1  and  2  (Fig  13B).  These  mutants  will  be  tested  for  potential  defects  in 
oligodendrocyte  and  cardiovascular  development,  along  with  the  mutant  lines  from  a  reverse  genetics 
approach  using  zinc-finger  nucleases,  discussed  below. 

Zinc  Finger  Nuclease  (ZFN)-induced  nfla  and  nflb  mutant  alleles 
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As  previously  outlined  in  our  proposal,  the  generation  of  stable  nfla  and  nflb  mutant  zebrafish 
lines  remains  a  necessary  prerequisite  to  performing  the  large  scale  chemical  and  genetic  screens  we 
argue  will  be  invaluable  to  furthering  our  understanding  of  NF1  biology  and  lead  to  the  development  of 
novel  therapeutic  agents.  While  we  encountered  some  success  in  generating  nfla  mutant  lines  using  the 
TILLING  approach,  this  method  lacked  the  efficiency  necessary  to  generate  multiple  loss-of-function 
mutant  alleles.  Therefore,  we  turned  to  emerging  technologies  that  promised  a  more  efficient  rate  of 
specific  lesion-induction  in  the  zebrafish  genome. 

Excitingly,  a  cutting-edge  strategy  for  targeted  mutagenesis  was  developed  for  zebrafish  using 
zinc  finger  nuclease  (ZFN)  technology.  We  were  fortunate  to  have  the  opportunity  to  collaborate 
closely  with  the  laboratories  of  Nathan  Lawson  and  Scot  Wolfe  to  apply  this  technology  to  generate 
stable  mutant  nfla  and  nflb  zebrafish  lines.  Briefly,  we  utilized  a  modular  approach  to  generate  two 
zinc  finger  proteins,  each  with  a  9-bp  sequence  specific  recognition  element.  Each  of  these  proteins  is 
coupled  to  a  monomer  of  the  heterodimeric  FokI  restriction  endonuclease.  In  vitro  synthesized  RNA 
encoding  these  ZFNs  are  injected  into  1-cell  zebrafish  embryos  where  they  are  translated  into  protein 
and  bind  to  there  cognate  target  sequences  in  the  zebrafish  genome  within  exons  of  nfla  and  nflb.  This 
binding  occurs  such  that  the  two  FokI  restriction  endonuclease  monomers  are  brought  into  close 
proximity  of  each  other  allowing  them  to  form  a  functional  dimer  and  introduce  a  double-stranded  break 
in  the  DNA.  These  double-stranded  breaks  are  repaired  by  the  error-prone  non-homolgous  end-joining 
(NHEJ)  pathway  leading  to  microinsertions  and  microdeletions  at  the  site  of  the  double-stranded  break. 
The  Look  and  Epstein  laboratories  have  collaboratively  generated  multiple  ZFNs  directed  against  nfla 
(exon  25)  and  nflb  (exon  19).  This  combined  effort  has  led  to  the  generation  7  different  nfla  and  nflb 
mutant  alleles  all  harboring  frame  shifts  that  result  in  premature  stop  codons  (Figure  14).  As  the  GAP- 
related  domain  in  nfla  and  nflb  span  exon  26  and  exon  27,  we  expect  that  these  frame-shift  mutations 
within  exon  25  of  nfla  and  exon  19  of  nflb  should  result  in  proteins  that  lack  functional  GAP  activity. 

Redundant  function  of  the  nfla  and  nflb  alleles 

We  generated  single  homozygous  mutant  alleles  of  nfla  and  nflb,  however  they  failed  to  show 
any  significant  abnormalities  from  embryogenesis  through  adulthood.  However,  when  n/Y/> specific 
morpholinos  were  injected  into  nflaZFN-mut2/ZFN-mut2  mutant  embryos,  an  increased  OPC  phenotype, 
similar  to  nfla  and  nflb  double  morphants,  was  observed  (data  not  shown).  This  result  suggests  that 
nflaZFN-mut2  has  no  functional  nfla  protein  and  that  nflb  shares  redundant  function  with  nfla. 

Due  to  the  lack  of  a  mutant  phenotype  nfla  and  nflb  single  homozygous  mutants,  we  generated 
nfla; nflb  double  mutants.  We  crossed  two  ZFN-induced  nfla  mutant  alleles  and  the  nfla- nonsense 
tilling  mutant  with  two  mutant  alleles  of  nflb  (nflaZFN-mut2;nflbZFN~mut2,  nflaZFN~mutl  ;nflbZFN~mut3  and 
nfla7203 ;nflbZFN~mut2).  To  analyze  tissue-specific  phenotypes  we  intercrossed  6  different  transgenic  lines 
with  the  nfla;nflb  double  mutants  (see  Table  2)  to  ascertain  cell  type-specific  phenotypes  in  the  CNS 
and  PNS.  We  also  found  that  nfla; nflb  double  homozygous  mutants  start  to  die  around  7  dpf  and  that 
the  majorty  fail  to  survive  beyond  9  to  10  dpf.  The  double  mutants  exhibit  severe  brain  degeneration  that 
may  be  a  primary  cause  of  lethality.  All  three  combinations  of  nfla;nflb  double  mutant  alleles  exhibited 
the  identical  phenotype  supporting  the  redundant  roles  of  nfla  and  nflb  and  that  they  are  essential  for 
survival  as  in  mammals. 

Hyperplastic  OPCs  and  SCG  in  the  nfla;nflb  double  mutant 

Although  nfla; nflb  double  homozygous  mutant  larvae  do  not  survive  beyond  10  dpf,  their 
general  morphology  appears  relatively  normal  until  7  or  8  dpf  (see  Fig.  15B  of  confirmed  mutant 
genotypes).  Most  double  mutant  larvae  lacked  melanocytes  in  the  lateral  stripe,  which  could  first  be 
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detected  at  4  dpf  and  persisted  until  death  (Fig.  15C  and  15D).  Furthermore,  the  double  mutant  larvae 
exhibited  increased  yellow  pigmentation,  indicative  of  altered  xanthophore  development  and  may 
correlate  with  increased  cafe-au-lait  spots  often  seen  in  human  NF1  patients.  We  have  now  confirmed 
that  we  can  use  this  phenotype  to  identify  double  mutant  animals  prior  to  PCR  genotyping,  which 
simplifies  our  experimental  protocols  immensely. 

Because  nfla  and  nflb  double  morphants  have  excess  OPCs  in  the  spinal  cord,  we  promptly 
confirmed  this  phenotype  in  Tg(olig4:egfp);nflaZFN-mut2ZZFN-mut2  ;nflbz mutants  and  found 
that  they  had  more  dorsally  migrated  OPCs  in  the  spinal  cord  than  wild-type  control  animals  at  6  dpf 
(Fig.  4E  and  4F).  Increased  numbers  of  OPCs  could  be  detected  as  early  as  3  dpf  in  the  double  mutants, 
which  persisted  up  to  their  death  (data  not  shown).  This  result  is  consistent  with  our  morphant  data  and 
the  mouse  Nfl-loss  phenotype. 

The  earliest  report  of  Nfl  knockout  mice  demonstrated  the  Nfl  homozygous  mutant  had 
hyperplasia  in  the  neurons  of  sympathetic  ganglion  tissue.  To  check  whether  zebrafish  nfla  and  nflb- 
loss  induces  a  similar  defect,  we  examined  superior  cervical  sympathetic  ganglion  (SCG)  tissue  by 
crossing  the  mutant  to  the  Tg(dbh:gfp)  transgenic  line.  As  shown  in  Fig  15G  and  15H,  the 

nflaZFN_mut2/ZFN_mut2;nflbZFN_mut2/ZFN_mut2  ^  an  enlaj.ged  SCG  relative  to  wild  type  animals, 

indicating  that  the  function  of  nfl  genes  in  zebrafish  is  conserved  with  Nfl  function  in  mice. 

In  order  to  address  whether  Ras  signaling  is  activated  in  the  double  mutant,  we  performed 
immunohistochemistry  using  phospho-ERK  1/4  antibody.  At  5  dpf,  a  few  neurons  in  the  spinal  cord  of 
wild-type  larvae  express  pERKl/4  strongly,  while  most  other  cells  express  it  weakly  (Fig.  151). 
However,  in  the  double  mutant,  in  addition  to  a  subset  of  neurons,  the  white  matter  of  neurons  in  the 
spinal  cord  are  show  strong  expression,  indicating  that  ERK  signaling  is  highly  active  in  the  axons  of 
neurons  in  the  mutant  (Fig.  15J).  Furthermore,  the  basal  level  of  ERK1/4  phosphorylation  appears 
relatively  higher  in  most  tissues  (Fig.  15J).  These  findings  suggest  that  nfla  and  nflb  loss  leads  to 
activation  of  the  ERK  pathway,  which  is  one  of  the  main  downstream  signaling  of  Ras  pathway. 

In  summary 

•  Transient  loss-of-function  of  nfla  and  nflb  in  zebrafish  yields  a  range  of  cardiac,  vascular,  and  neural 
phenotypes  that  we  have  carefully  defined  and  characterized. 

•  Several  stable  lines  of  zebrafish  harboring  loss-of-function  alleles  for  nfla  and  nflb  have  been 
identified  and  validated 

Materials  and  Methods 

1.  OPC  Imaging  and  statistical  analyses.  For  time-lapse  movies,  a  Zeiss  200M  inverted 
microscope  equipped  with  a  spinning-disc  confocal  system  (Yokogawa)  was  used.  The  statistical 
significance  of  changes  to  OPC  cell  numbers  and  migration  upon  nfl  knockdown  was  determined  by  t- 
test  analysis. 

2.  GRD  rescue.  The  construct  for  containing  the  mCherryRed-GRD  fusion  protein  was  first 
created  using  the  multisite  Gateway  system  (InVitrogen)  and  subcloned  into  the  pCS2+  vector  in  order 
to  synthesize  sense  RNA  in  vitro.  Sense  RNA  for  the  mCherryRed  fluorescent  protein  was  synthesized 
and  served  as  an  injection  control.  Approximately  15  pg  of  mCherryRed-GRD  RNA  or  6  pg  of 
mCherryRed  were  injected  for  the  rescue  experiments,  along  with  appropriate  morpholinos. 

3.  In  situ  hybridization,  nfla  (ENSDARG000000 12982)  and  nflb  (ENSDARG00000004184) 
predicted  cDNA  sequences  from  Ensembl  (http://www.ensembl.org/)  were  subjected  to  VISTA  analysis 
(http://genome.lbl.gov/vista/)  and  the  two  regions  of  greatest  dissimilarity  were  PCR  amplified  for  both 
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genes  (nfla  from  bp  865-1518:  5 '-GGGGG  AG  AGG  AG  A  A  AGCT  A  A  and  5'- 

GGGGTCAGCGTGGATAAGT ;  nfla  from  bp  6638  to  7236:  5'-CCACAGTCATCGCTCTCACC  and 
5 '-GGCT ACTGCGTGAAAACC AG;  nflb  from  bp  770  tol491:  5 '-CCGAAAGC ACGAAG AGAAAA 
and  5 '-TCTGCTCTTGGCTCTGGTG;  nflb  from  bp  6798  to  7320:  5'- 

CGAGGAGGTGTTTATGGAGA  and  5  '-GGTGTTGACTTGAGAGGGTTG)  from  24  hpf  pooled 
zebrafish  cDNA  and  cloned  into  the  pCRII-TOPO  vector  (Invitrogen).  GFAP  in  situ  probe  was  PCR 
amplified  from  zebrafish  3dpf  pooled  cDNA  (5'-  AGCGTTCCTTCTC ATCCT ACCGAA  and  5'- 
GCAGTCTCGTCTTGAAGTCTCTGT)  and  cloned  into  the  pCRII-TOPO  vector  (Invitrogen). 

Antisense  RNA  probes  for  nfla  and  nflb  (from  bp  6798-7320)  were  generated  by  linearization  with 
EcoRV  digestion,  followed  by  transcription  with  SP6  polymerase.  An  antisense  RNA  probe  for  nflb 
(from  bp  770  to  1491)  was  generated  by  linearization  with  Notl  digestion,  followed  by  transcription  with 
SP6  polymerase.  Processed  embryos  were  mounted  in  100%  glycerol  and  photographed. 

4.  Western  blotting.  Whole  zebrafish  3.5  dpf  protein  lysates  (50-100pg)  were  separated  by  gel 
electrophoresis,  transferred  to  PVDF  membranes  (Invitrogen)  and  probed  overnight  at  4°C  with  the 
following  primary  antibodies:  anti-neurofibromin  (Abeam  Abl7963,  1:1000),  anti-phospho-p44/42 
MAPK  (Cell  Signaling  7100,  1:100;  Cell  Signaling  4377,  1:1000),  anti-p44/42  MAPK  (Cell  Signaling 
9102,  1:1000),  and  anti-alpha-tubulin  (Sigma  T6074,  1:1000).  Primary  antibody  binding  was  visualized 
using  the  Western  Breeze  Immunodetection  Kit  (Invitrogen)  according  to  manufacturer  instructions. 

5.  Imaging  and  statistical  analysis.  Developing  Tg(flil:egfp)yl,Tg(fli:negfp)y7,  and 
Tg(kdr:GRCFP)znl  zebrafish  embryos  were  mounted  in  low-melt  agarose  supplemented  with  tricaine 
methanesulfonate  (MS-222)  and  imaged  using  a  Leica  TCS  SPII  confocal  laser  scanning  microscope. 
Single  vertical  projections  of  three-dimensional  datasets  (60  optical  sections  at  2pm  intervals  from  a 
10X  water  immersion  lens)  were  generated  as  TIFF  images  using  the  Leica  Confocal  Software  package 
(Leica).  Transmitted  light  images  were  captures  using  an  Olympus  MVX10  microscope  and  Olympus 
DP72  digital  camera  (Olympus).  The  SAS/STAT  software  package  (SAS)  was  used  to  perform 
statistical  analyses  for  ordinal  data  using  a  Mann- Whitney  U  test  while  continuous  data  was  analyzed 
using  a  Student's  t-test. 

6.  Histology.  Briefly,  samples  were  washed  with  ice-cold  PBS,  fixed  overnight  in  4% 
paraformaldehyde,  washed  with  PBS,  and  dehydrated  through  an  ethanol  series  prior  to  paraffin 
embedding.  Anti-phospho-p44/42  MAPK  antibody  (Cell  Signaling  4370)  was  used  at  manufacturer 
suggested  concentration. 

7.  Flow  Cytometry  and  Quantitative  PCR.  Fluorescence  activated  cell  sorting  (FACS)  was 
performed  on  whole  single  cell  suspensions  generated  from  3.5  dpf  wild  type  and  Tg(kdrl:GRCFP)znl 
zebrafish  embryos.  Briefly,  embryos  were  rinsed  in  calcium  free  Ringer’s  solution  (116  mM  NaCl,  2.9 
mM  KC1,  5  mM  HEPES  [pH=7.2]),  deyolked,  and  incubated  in  protease  solution  (0.25%  trypsin,  1  mM 
EDTA;  in  PBS  [pH=8])  at  28  degrees  Celsius  with  intermittent  agitation  by  repeated  pipetting  through  a 
200  uL  tip.  Upon  generation  of  a  single  cell  suspension,  reaction  was  stopped  by  addition  of  10%  fetal 
calf  serum.  Cells  were  rinsed  in  suspension  solution  (Leibovitz  L-15  Medium,  0.8  mM  CaC12,  Penicillin 
[50  U/mL],  Streptomycin  [0.05  mg/ml],  1%  fetal  calf  serum),  filtered  successively  through  100  pm  and 
40  pm  cell  strainers,  and  sorted  using  a  BD  FACS  Aria  cell  sorter  into  collection  medium  (Leibovitz  L- 
15  Medium,  0.8  mM  CaC12,  Penicillin  [50  U/mL],  Streptomycin  [0.05  mg/ml],  10%  fetal  calf  serum). 
Quantitative  PCR  was  performed  for  cdh5  (5'-AAATGGCACAGACCTCTTCA  and  5'- 
CTCACGGTCGAGGTCACTT),  etsrp  (5 ACTCCTGCTGG ATTCTGCTT  and  5'- 
GGCCACGACTGAGCTTCT),  and  efla  (5 '-CTTCGTCCC AATTTC AGGAT  and  5'- 

GAACC  AGCCC  ATGTTTGAG) . 
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Key  research  accomplishments: 

•  Characterization  of  zebrafish  orthologs:  nfla  and  nflb. 

•  Further  Analysis  of  gene  expression  ofnfl  a  and  nfl  b. 

•  Further  analysis  ofmorpholino  knockdown  phenotypes  for  nfl  a  and  nfl  b. 

•  Generation  of  zebrafish  lines  harboring  mutations  in  nfla  and  nflb  via  TILLING  and  through 
the  use  of  engineered  Zinc  Finger  Nucleases. 

•  Initial  phenotypic  analysis  ofnfla/nflb  double  homozygous  mutant  zebrafish  embryos. 

Reportable  outcomes: 

A  manuscript  reporting  our  cardiovascular  findings  following  MO-mediated  transient  loss-of-function  of 
nfla  and  nflb  was  published  in  PNAS. 

Cardiac  and  vascular  functions  of  the  zebrafish  orthologues  of  the  type  I  neurofibromatosis  gene  NFI. 
Padmanabhan  A*,  Lee  JS*,  Ismat  FA,  Lu  MM,  Lawson  ND,  Kanki  JP,  Look  ATA,  Epstein  JAA. 

Proc  Natl  Acad  Sci  USA.  2009  Dec  29;106(52):22305-10.  Epub  2009  Dec  4. 

*  =  These  authors  contributed  equally  to  this  work 
A  =  Corresponding  authors 

A  second  manuscript  reporting  our  neural  findings  following  MO-mediated  transient  loss-of-function  of 
nfla  and  nflb  is  presently  in  the  review  process  at  Human  Molecular  Genetics.  In  addition,  multiple 
abstracts  and  oral  presentations  stemming  from  this  work  have  been  presented  at  various  meetings 
including  the  University  of  Pennsylvania  Cardiovascular  Institute  Symposium,  University  of 
Pennsylvania  Cell  and  Molecular  Biology  Symposium,  The  6th  European  Zebrafish  Development  and 
Genetics  Meeting,  and  The  Mid- Atlantic  Regional  Zebrafish  Meeting. 

Conclusion: 

In  this  third  year  of  study,  we  have  significantly  furthered  our  advances  towards  meeting  the  goals  set 
forth  in  the  funded  proposal.  The  target  genes  have  been  clearly  identified  and  phenotypes  due  to 
morpholino  knockdown  have  been  characterized  that  are  relevant  to  the  understanding  of  human  disease. 
Our  efforts  towards  this  end  have  resulted  in  a  published  manuscript  with  a  second  presently  in  the 
review  process.  While  morpholino  knockdown  is  informative,  future  experiments  will  be  aided  by  the 
characterization  of  stable  mutations  within  these  genes  in  zebrafish.  We  have  now  obtained  these 
mutant  fish  lines  using  several  different  approaches.  When  characterized,  these  reagents  will  be  freely 
shared  with  our  colleagues.  These  reagents  should  be  excellent  candidates  for  use  in  screening  small 
molecules  and  for  use  in  performing  genetic  screens  to  identify  modulators  of  this  important  disease. 


Page  1 1 


Supporting  Data 


Chr  15 
Forward  strand 

I - H-H-HHH 

nfla  genomic 


-m 


10 

iiiiiiiiiiiiiiiai 

nfla  mRNA 


20 

iiiii 


nfla 

240  kB 
120  kB 

— H-t~ 


I  • 


\\\h 


30 

III  lllllll  IIIII 

Ras  GRD 
8268  bp 


B 


40  50 

lll■llllllllllllllll 


Chr  10 
Forward  strand 


nflb 

270  kB 
134  kB 


4HHI  IIH  IH  I  H  IM 


nflb  genomic 


IH  HH-hl-HH-fl 


1  10  20  30 

llllllllllllllll■ll■lllll  lllllll  II 

nflb  mRNA  Ras  GRD 

- 8244  bp 


40  50 

llll■■lllllllllllll■l 


-l' 


-  drosophila  neurofibromin 


bOrafir-lk  neurof'broin*!*  b 
human  neurofibromin 
mouse  neurofibromin 
zebrafish  neurofibromin  a 


*  Is 

5  3* 


*  *  s 


5  8 


MS/i’ 


hChr  17 


msi2  I  zChr  10 
25  Kb 


o 


Fig  1.  Zebrafish  have  two  orthologues  of  human  NF1. 

(A,  B)  Genomic  and  mRNA  structures  of  the  two  orthologous  zebrafish  genes  corresponding  to  human 
NFL  (C)  Phylogenetic  tree  comparison  of  zebrafish,  human,  mouse,  and  Drosophila  neurofibromin.  (D) 
Analysis  of  syntenic  relationships  between  human  chromosome  17  (NF1)  and  zebrafish  chromosomes 
15  (nfla)  and  10  (nflb).  Relative  genomic  positions  are  to  scale  as  indicated. 
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ETVLADRFERLVELVTMMGDQGELPIAMALANVVPC50WDELARVLVTLFDSRHLLYQLLWNMFSKEVELADSMQTLFRGNSLASKIMTF 

etvladrferlvelvtmmgdqgelpiamalanvvpcbqwdelarvlvtlfdsrhllyqllwnmfbkeveladsmqtlfrgnslaskimtf 

ETVLADRFERL VELVTMMGDQGELP I AMALANV VPC6QWDELARVLVTLFDSRHLLYQLLWNMFBKEVELAD5MQTLFRGNSLASK IMTF 

etaladrferlvelvtmmgdqgelpiamalasvvpcbqwdelarvlvtlfdsrhllyqllwnmfskbveladsmqtlfrgnslaskimtf 

etvladrfeqlvqlvtmisdkgelpiamalanwttsqmdelarvlvtlfdakhllspllwnkfyrevevsdcmqtlfrgnslgskimaf 

CFKVYGATYLQKLLDPLLRI V I T S S DWQH VS FE VDPTR LEP SE S LEE NQRN LLQMTE K F F H A I I SS S S E F P PQLRS VC HC L Y Q ATC H S LL 
cfrvygatylqklldpllrvi ITSSDWQHVSFEVDPTP.LEPSESLEENQRNLLQMTERFFHAI issssefpsqlrsvchclyqatchsll 
CPKVYGAAYLQKLLEPLLRGI ITTPEWQQXSFE VDPTFLBQAENLEENQKNLLQITDRFFLAI IhSSSEFPPQLRSVCHCLYQATCHSLL 
CFKVYGATYLQKLLEPLLRSVIAGTEWHMLSFEVDPTRLENGENLEENQRNLLRMTERFFKAIISSSGEFPPQLRSVCHCLFQATCHSLL 
CFKIYGASY  LQMLLE PLIR PLLD--EE E ETCFE VDP AR LD PTED I EQHRNNL.  I ALTQKVF DAI INSSDRFPPQLRSMCHCLYQ - 

MKATVXEKKENKKSVVSQRFP QNS IG AVGS AM FLR F I N P A I VSP YE AG1 L DKK P P PR  I E RGLK LM S K I LQS I  A N H VLPTK E E H MR P F 

NKATVKERKENKKSVVSQRFP QNS IGAVGS AMFLR  FIMPAI VSPYE AGILDKKPPPR IERGLKLMSKVLQSIANHVLPTKEEHMRPF 

NKATVKEKKENKKSVVSORFP QNS IGAVGS AMFLR FI N P A I VS P YE AGILDKK PPPR I E RGLK LMS K I LQS I ANH VLFT KEEHMR P F 

SKASVKDRKEVRKAVVSQRFP ONS I G AVGS AM F LR FVN P A I VS P YE AG I LD K K PQPR I E RGLK LM S K I LQS I  A N H VLFT KE E HMR P F 

- VLSKRFPNLLQNNIGAVGTVI F LRFI N P A I VSPQELGI VDKQ VHS S AK RGLMLMSK I LQN I ANHVE FSKEQHMLCF 

HDFVKSNFDAARRFFLDIASDCPTSDAVNHSLSFISDGNVLALHRLLWNNOEKIGQYLSSNRDHKAVGRRPFDKMATLLAYLGPPEH 

NDFVKSNFDLARRFFLDIASDCPTSDAVNHSLSF1SDGNVLALHRLLWNN0EKIGGYLSSNRDHKAVGRRPFDKMATLLAYLGPPEH 

NDFVKSNFDAARRFFLD1ASDSPPBDSVNHSLSFISDGNVLALHRLLWNN0ER1G0YLS5NRDHKAVGRRPFDKMATLLAYLGPPEH 

ndfvksnfdaarrffldiasdcpabdslnhslsfisdgnvlalhrllwnnqekigqylsbnrdhkavgrrpfdkmatllaylgppeh 

NDFLRDHFEAGRRFFIGIASDCETV00T3HSMSFISDANVLALHRLLWTHQEKI G  D YLS  SSRDHK A VGRRP  FDKMATLL A  Y  LG  P  PEH 


Fig  2.  Zebrafish  NF1  orthologues  are  highly  conserved  within  the  GAP-related  domain. 

Amino  acid  sequence  alignment  reveals  a  remarkably  high  degree  of  protein  conservation  within  the 
human,  mouse,  zebrafish,  and  drosophila  neurofibromin  GAP-related  domains  (identical  amino  acid 
highlighted  grey,  similar  amino  acids  highlighted  yellow,  and  mismatched  amino  acids  not  highlighted). 
Alternatively  spliced  exon  23a  in  human  NF1  (not  depicted)  is  present  in  both  nfla  and  nflb. 
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Fig  3.  nfla  and  nflb  are  expressed  maternally  and  in  the  developing  zebrafish  cardiovascular 
system. 

Whole-mount  in  situ  hybridization  for  nfla  and  nflb  at  the  4-cell  stage,  24  hpf,  48  hpf,  3dpf,  and  4  dpf. 
(Al,  Bl)  At  the  4-cell  stage,  nfla  and  nflb  are  expressed  throughout  the  animal  pole  of  the  developing 
embryo.  (A2,  B2)  Both  genes  are  expressed  broadly  at  24  hpf  (inset),  with  strong  expression  along  the 
spinal  cord.  (A3,  B3)  At  48  hpf,  expression  of  nfla  and  nflb  is  noted  in  the  head  and  regions  of  the 
anterior  trunk  (inset).  Spinal  cord  expression  of  both  genes  persists  and  positive  staining  is  observed 
along  the  dorsal  vessel  for  nfla  and  nflb.  (A4,  B4)  Cardiac  expression  for  both  genes  is  observed  at  48 
hpf.  (A5,  B5)  Expression  of  nfla  and  nflb  become  progressively  restricted  to  regions  of  the  head  at  3 
dpf  (insets),  nfla  and  nflb  expression  along  the  dorsal  vessel  (A5,  B5)  and  in  the  embryonic  heart  (A6, 
B6)  persist  at  3  dpf.  (A7,  B7).  At  4  dpf,  robust  vascular  staining  is  apparent  for  nfla  and  nflb.  Scale 
bars:  25pm  (100pm  for  insets). 
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Fig  4.  Zebrafish  nfla  and  nflb  are  expressed  during  embryonic  development. 

(A)  Distribution  of  ESTs  mapping  to  nfla  or  nflb  in  dbEST  including  tissue  of  origin.  (B)  RT-PCR 
analysis  verifying  maternal  (1-cell)  and  zygotic  (24,  72,  and  84  hpf)  expression  of  nfla  and  nflb.  (Cl- 
C6,  D1-D6)  Whole-mount  in  situ  hybridization  of  24  hpf,  48  hpf,  3  dpf  (all  magnifications  of  insets  in 
Fig.  3),  and  4  dpf  zebrafish  embryos  shows  temporal  and  spatial  expression  patterns  of  nfla  and  nflb. 
(C5-6,  D5-6)  Both  nfla  and  nflb  are  expressed  in  cardiac  tissue  at  48  hpf  and  3  dpf  as  seen  in  these 
anteroventral  images.  Notably,  analysis  of  14pm  (C7,  D7)  cross  sections  of  whole-mount  stained  4  dpf 
zebrafish  embryos  reveals  strong  nfla  and  nflb  expression  in  the  region  of  the  dorsal  vessel  (arrows) 
directly  below  the  notochord  (NC).  (E)  Parameters  utilized  for  sorting  GFP-positive  and  GFP-negative 
cells  from  3.5  dpf  Tg(kdrl:GRCFP)znl  zebrafish  embryos.  (F)  Quantitative  PCR  analysis  of  the  GFP- 
positive  cell  fraction  reveals  marked  enrichment  for  endothelial-specific  genes,  including  cdh5  and 
etsrp,  relative  to  the  GFP-negative  fraction.  (G)  Reverse-transcription  and  polymerase  chain  reaction 
(RT-PCR)  for  nfla  and  nflb  was  performed  using  first  strand  cDNA  synthesized  from  endothelial  cell 
RNA  (GFP-positive  cell  fraction),  confirming  expression  of  both  genes  in  vascular  endothelium.  Scale 
bars:  25pm. 
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Fig  5.  Histological  analysis  of  MO  knockdown  cardiac  valves. 

Serial  transverse  histological  sections  through  the  atrio-ventricular  valves  of  wild  type  (Al-4)  or  ~2ng 
morphant  nfla  ATG  (Bl-4),  nflb  ATG  (Cl-4),  and  nfla  +  nflb  ATG  (Dl-4)  3.5  dpf  zebrafish  embryos 
reveal  no  readily  apparent  defects  at  the  resolution  afforded  to  us  by  histological  analysis  (A,  atrium;  V, 
ventricle).  Scale  bars:  50pm. 
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Fig  6.  MO  knockdown  of  nfla,  nflb,  or  both  together  results  in  pericardial  effusions  at  3.5  dpf  and 
increased  phospho-p44/42  MAPK  in  cardiac  tissue. 

Analysis  of  3.5  dpf  wild  type  (A),  nfla  ATG  5-mispair  (5MP)  MO-  (B),  nflb  ATG  5MP  MO  (D),  or 
nfla  +  nflb  ATG  5MP  MO-  (F)  injected  embryos  reveal  no  apparent  defects  in  gross  morphology. 
Treatment  with  nfla  ATG  MO  (C),  nflb  ATG  MO  (E),  or  a  combination  of  both  (G),  however,  results 
in  a  dilation  of  the  pericardial  space.  (H-J)  Injection  of  p53^  embryos  with  nfla  +  nflb  ATG  MO  results 
in  a  gross  dilation  of  the  pericardial  space  (J),  while  uninjected  (H)  and  nfla  +  nflb  ATG  5MP  MO- 
injected  p53~'~  embryos  (I)  appear  normal.  Scale  bars:  0.25mm.  (K-L)  Transverse  sections  of  3.5  dpf 
nfla/nflb  combined  morphant  embryos  reveals  a  thinning  of  the  ventricular  myocardium  and  pericardial 
effusion  (*)  when  compared  with  controls  (A,  atrium;  V,  ventricle).  (M-N)  Immunohistochemical 
analysis  of  transverse  sections  of  3.5  dpf  nfla/nflb  combined  morphant  embryos  reveals  an  increase  in 
the  ratio  of  phospho-p44/42  MAPK  positive  cardiac  cells  (arrows)  to  the  total  number  of  cardiac  cells 
when  compared  to  controls.  Scale  bars:  25pm. 
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Fig  7.  MO  knockdown  of  nfla  or  nflb  results  in  a  thinning  of  the  ventricular  myocardium  and 
increased  phospho-Erk  staining  in  cardiac  tissue. 

(A,  A',  B,  B',  F,  G)  Transverse  histological  sections  through  the  heart  of  3.5  dpf  zebrafish  embryos 
reveals  a  thinning  of  the  ventricular  myocardium  (red  arrows)  of  zebrafish  treated  with  nfla  or  nflb 
ATG  MOs  when  compared  with  the  corresponding  controls  (A,  atrium;  V,  ventricle;  YS,  yolk  sac). 
Notably,  a  large  pericardial  effusion  can  be  appreciated  in  nfla  and  nflb  morphants  (*).  (C,  D,  H,  I) 
Immunohistochemichal  analysis  of  transverse  sections  through  3.5  dpf  zebrafish  hearts  reveals  an 
upregulation  of  phospho-p44/42  MAPK  (phospho-Erk)  in  nfla  or  nflb  morphants  when  compared  to  the 
appropriate  5MP  controls.  The  ratio  of  phospho-Erk  positive  cardiac  cells  to  the  total  number  of  cardiac 
cells  was  quantified  over  multiple  sections  for  nfla  MO  (2ng),  nflb  MO  (2ng),  and  nfla  +  nflb  MO 
(2ng)  groups  (mean  ratio  phospho-p44/42  MAPK  positive  cardiac  cells/total  cardiac  cells  +  SD)  (E,  J, 
K).  (L)  Quantification  of  percentage  of  embryos  displaying  a  pericardial  effusion  phenotype  at  3.5  dpf 
following  injection  of  approximately  2ng  of  nfla  ATG  5MP  MO,  nfla  ATG  MO,  nflb  ATG  5MP  MO, 
or  nflb  ATG  MO.  75%  of  nfla  ATG  morphants  displayed  pericardial  effusions  as  compared  to  5%  of 
nfla  ATG  5-mispair  (5MP)  morphant  controls  (p<0.0001).  24%  of  nflb  ATG  morphants  displayed 
pericardial  effusions  as  compared  to  4%  of  nflb  ATG  5MP  morphant  controls  (p<0.01).  (M) 
Quantification  of  percentage  of  embryos  displaying  a  pericardial  effusion  phenotype  at  3.5  dpf  following 
injection  of  nfla  ATG  MO,  nflb  ATG  MO,  or  nfla  +  nflb  ATG  MO.  A  suboptimal  dose  (lng)  of  nfla 
ATG  MO  resulted  in  a  12%  incidence  of  pericardial  effusion,  while  a  suboptimal  dose  (lng)  of  nflb 
ATG  MO  resulted  in  a  7%  incidence  of  pericardial  effusion.  Injection  of  lng  of  a  combination  of  nfla  + 
nflb  ATG  MO,  however,  resulted  in  a  20%  incidence  of  pericardial  effusion.  Scale  bars:  25pm. 


Page  18 


Fig  8.  MO  knockdown  of  nfla  results  in  vascular  patterning  defects  at  48  and  72  hpf. 

(A-C)  At  48  hpf,  nfla  ATG  MO-treated  Tg(fli:egfp)yl  (endothelial-specific  GFP  transgenic)  zebrafish 
embryos  display  gross  defects  in  vascular  development  compared  with  control  MO-treated  or  uninjected 
samples.  Morphant  embryos  (C)  display  abnormal  claw-like  projections  at  the  leading  edge  of  the 
developing  intersomitic  vessels  and  fail  to  develop  the  dorsal  longitudinal  anastomotic  vessel  (DLAV) 
present  in  both  control  MO-treated  (B)  and  uninjected  (A)  samples.  (D-F)  At  72  hpf,  nfla  ATG 
morphant  embryos  display  only  rudimentary  DLAVs  and  a  general  disorganization  of  the  trunk 
vasculature  (F)  when  compared  with  control  MO-treated  (E)  or  uninjected  (D)  embryos.  Scale  bars: 

25  pm. 
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Fig  9.  MO  knockdown  of  nfla,  nflb,  or  both  together  result  in  vascular  defects  at  24  at  48  hpf. 

(A,  B)  Analysis  and  quantification  of  vascular  defects  at  24  hpf  in  uninjected  and  morphant 
Tg(fli:negfp)y7  (endothelial-specific  nuclear  GFP  transgenic)  zebrafish  embryos.  Control  MO-  (Al)  and 
combined  nfla/nflb  MO-treated  (A2)  zebrafish  embryos  appear  similar  by  gross  morphological  analysis 
at  24  hpf  (scale  bars:  500pm).  Development  of  intersomitic  vessels  is  deficient  at  24  hpf  in  nfla/nflb 
combined  morphants  (A4)  when  compared  to  controls  (A3)  (scale  bars:  25pm).  (B)  Intersomitic  vessel 
formation  between  somites  17-30  at  24  hpf  was  scored  as  absent  (red),  intermediate  (grey),  or  normal 
following  administration  of  2ng  of  the  indicated  MO(s).  (C)  MO-mediated  knockdown  of  flt4,  providing 
a  sensitized  background  for  the  detection  of  vascular  defects,  was  combined  with  nfla,  nflb,  and  nfla  + 
nflb  ATG  MO  knockdown.  24%-85%  of  combined  flt4/(nfl a ,  nflb,  nfla  +  nflb )  MO-treated  embryos 
display  abnormal  vascular  shunts  compared  with  3-8%  of  individual  flt4,  nfla,  nflb,  or  nfla  +  nflb  MO- 
treated  embryos. 
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Fig  10.  MO  knockdown  of  nfla,  nflb,  or  both  together  result  in  vascular  defects  at  24  at  48  hpf. 

(A)  Quantification  of  vasculature  at  24  hpf  in  uninjected  and  morphant  Tg(fli:negfp)y7  (endothelial- 
specific  nuclear  GFP  transgenic)  zebrafish  embryos.  Uninjected,  nfla- SBel  5MP  (lng),  and  nfla- SBel 
(lng)  morphant  embryos  were  qualitatively  assessed  for  presence,  absence  (red),  or  an  intermediate 
phenotype  (grey)  with  regard  to  the  developing  trunk  intersomitic  vessels  between  somites  17-30  at  24 
hpf.  (B)  Quantitation  of  vascular  shunt  phenotype  seen  with  low  (lng)  doses  of flt4,  nfla,  nflb,  or  nfla 
+  nflb  MOs.  (C)  Quantitation  of  cystic  expansion  phenotype  observed  in  morphant  Tg(fli:negfp)y7 
embryos  treated  with  2ng  of  the  indicated  MO(s).  (D1-D8)  Treatment  of  Tg(fli:negfp)y7  zebrafish 
embryos  with  an  nfla  ATG  MO,  nflb  ATG  MO,  or  nfla  +  nflb  ATG  MOs  leads  to  a  cystic  expansion 
in  the  region  of  the  caudal  artery  and  vein  (arrows,  D2-D4)  not  present  in  uninjected  embryos  (Dl).  The 
expanded  tissue  was  confirmed  to  be  vascular  by  GFP  expression  (D6-D8).  Scale  bars:  50pm. 
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Fig  11.  nfl  loss  affects  the  migration  of  OPCs. 

(A,  C)  An  uninjected  olig2-EGFP;  p53  e7/e7  control  embryo;  (B,  D)  An  olig2-EGFP\  p53  e7/e7 
embryo  injected  with  nfla+lb  morpholinos.  (A,  B)  Montages  of  12-hour  time-lapse  images  from  60  hpf 
to  72  hpf,  showing  the  movements  of  OPCs  every  hour.  Numbers  in  each  panel  denote  the  hour(s)  after 
the  start  of  imaging.  For  the  complete  movies,  refer  to  Movie  SI,  S2.  (C,  D)  Cell  migration  paths  are 
shown  for  five  OPCs  that  traveled  the  farthest  during  the  observation  period.  In  both  conditions, 
migration  patterns  are  highly  dynamic.  Arrowheads  indicate  the  end  point  of  the  cell’s  migration.  The 
dashed  lines  represent  the  dorsal-most  GFP+  domains  of  the  ventral  spinal  cord.  Top,  dorsal  (d);  left, 
rostral  (r).  (E)  Graph  showing  the  total  distance  OPCs  travel  with  or  without  nfl  loss.  (F)  Graph 
showing  the  total  time  individual  OPCs  spent  pausing.  (G)  Graph  showing  the  velocity  of  OPCs 
calculated  from  when  they  were  actively  migrating.  All  individual  OPCs  (represented  by  single  points) 
that  could  be  observed  in  the  field  of  view  throughout  the  entire  period,  spanning  the  12-hour  imaging 
period,  were  traced  and  used  in  this  analysis  (n=9  for  uninjected  control  embryos,  n=23  for  nfla+lb 
morphant  embryos).  The  data  in  E-G  are  reported  as  mean  ±  SEM;  asterisks  indicate  nonparametric 
statistical  significance  (*  <  0.05;  ns,  not  significant).  Scale  bars  in  A,  B  =  25  pm,  C,  D,  the  X  axis  = 
10pm  and  the  Y  axis  =  5  pm. 
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Fig  12.  Forced  expression  of  the  GRD  of  NF1  is  sufficient  to  rescue  the  aberrant  phenotype. 

(A-D)  Projected  confocal  images  showing  the  lateral  view  of  the  spinal  cord  of  olig2-EGFP\  p53  e7/e7 
embryo  at  3dpf:  A,  B,  nfla+lb  morpholino;  C,  D,  control  morpholino;  A,  C,  GRD  RNA-injected;  B,  D, 
mCherryRed  (control)  RNA-injected.  Insets  show  the  bright  field  images  of  the  embryos  used  for 
confocal  imaging  and  demonstrate  normal  overall  morphologies  in  all  cases.  Scale  bar  =  50  pm.  (E) 
Graph  showing  the  number  of  OPCs,  represented  by  individual  points,  in  the  dorsal  spinal  cord  (SC)  of 
10  embryos  for  each  condition  at  3  dpf.  (*  <  0.05;  ***  <  0.001).  The  data  is  presented  as  means  ±  SEM. 
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CCG-^CTG 

Pro->Lys 

3  (dir. 
seq.) 

Exon26 

CAG->CAC 

Gln->His 

4(74-11) 

Exon26 

TCA-»CCA 

Ser-»Pro 

5  (VU180) 

Exon26 

ATG^AAG 

Met->Lys 
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Fig.  13.  Searching  for  zebrafish  mutants  for  nfla  and  nflb 

(A)  Schematic  of  TILLING  procedure.  (B)  A  list  of  the  mutations  for  nfla  and  nflb  identified  from 
TILLING  thus  far. 
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wt  nfla  <?xon25 
nfl£L_ZFN_mutl 
nfla_ZFN_mut2 


TTCCTCTACAGCCTCAGGAGGGT 

GTCTTCCT  CT  Af  AGCCTGAGGAGCC  T  GACGGT  GTTGAGCT  CAT  GGAGGCCAAAT  C  ACAGCTCTTCCTGAA 


GTcncc  icitc . - . -ctgaa 

GTCTTCCTCT - AGGAGGGTGACGGTCTTGAGCTCATGGAGGCCAAATCACAGCTCTTCCTCAA 


wt  nfla  exon25  I PLQPFFGDGVFl  MFAKSQl  FI 
nf la_ZFN_mutl  LPL<  ;  - 
nfla_ZFN  mut2  LPLGG 


wt  nflb  exon!9 
nflb__ZFNLmutl 
nflb_ZFN_mut2 
nfib_ZFN_mut3 
nflb_ZFNI_mut4 
nf lb_ZFN_mut5 


T  r CCG TT ACCT GT GT GAGCAGCCA 

TGGAGGCCGTTCTAGTGGCCATGTCCTGTnCCGTTACCTGTGTGAGGACGCAGAAATCCGGTGCGGTGTTGAAGACATCCCAGTGCAGTCTCTTaGCCAAA 
TGGAGGCCGTTCTAGTGGCCATGTCCTGTTTCCGT-  CC -TGT GAGGAGGC AGA AATCCGGT GCGGT GTT GAAGACAT CCCAGT GCAGT CT CTT CT GCC A A A 
1 GGAGGCCGT I  C1 AGIGGCCAT'GTCCTGr ( I CCGf I  ACC  Cl GAGGAGGC  AGAAA I CCGG I GCGG I G1 1 GAAGACAT CCCAG1 GCAGTCT Cl IC  IGCCAAA 

T  GGAGGCCGTT  CT  AGT  GGCCAT  GT  CCT  GTTT  CCCTTACC . - . . .  AA 

TCGAGGCCGTTCTAGTGGCCATGTCCTGTTTCCGTTA -  T  T  _  GAGGAGGC  AG AAAT CCGGT GCGGTGTTGA AG ACAT CCCAGT GCAGT CT CTT  CT GCCAAA 
T  GGAGGCCGTT  CT  AGT  GGCCAT  GTCCT  GTTT  C  C  GTT - CC - GT  T  . GAGGCA  GA AAT CCGGT GCGGT GTT  GAAGACAT CCCAGT GCAGT CT CTT CT GCCAAA 


wt  nflb  exonl9 
nflb_ZFN_mutl 
nflb_ZFN_mut2 
nflb_ZFN_mut3 
nflb_ZFN_mut4 
nflb_ZFN_mut5 


FAVl  VAMSCFRYl  CEFAFIRCGVFDTPVQSl  l  P 

EAVLVAMSCFRGL 

FAVLVAMSCFRYPFPFP' 

EAVLVAMSCFRYLPLTTTPSWSLPLSV7 

EAVLVAMSCFR 

EAVLVAM5CFRSVFRGRNPVRC* 


Fig.  14.  ZFN-induced  nfla  and  nflb  stable  mutant  zebrafish  lines 

Schematic  representation  of  the  nfla  and  nflb  mutant  alleles  generated  using  ZFN  technology. 
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nflaLtNmua; 

nflbZFNmut2 

nflaLtNmutl; 

nflbZFNmut3 

nflbZFNmut2 

marked  cell  types 

Tg(olig2:egfp) 

O 

O 

O 

OPCs  and  oligodendrocytes  (CNS) 

Tg(soxlO:egfp) 

O 

O 

X 

OPCs  and  oligodendrocytes 

(CNS),  Schwann  cells  (PNS) 

Tg(nkx2.2a:mgfp) 

O 

O 

X 

developing  myelin  (CNS), 
perineurial  glia  (PNS) 

Tg(  3. 9nestin:gfp ) 

O 

X 

X 

neural  precursor  (CNS) 

Tg(brn3c:gfp) 

O 

X 

X 

optic  nerve 

Tg(dbh:gfp) 

O 

X 

X 

sympathetic  ganglion 

Tg(kdrl:grcfp) 

O 

O 

X 

endothelial  cells 

Tg(fli:negfp ) 

O 

O 

X 

endothelial  cells  (nuclear) 

Tg(c236:GFP) 

X 

X 

X 

developing  AV  valve 

Table  1.  Three  nfla;nflb  double  mutant  alleles  crossed  with  transgenic  lines  (0=  line  is  generated;  X= 
line  is  being  established.) 
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Fig  15.  Phenotypes  of  nfla;nflb  double  mutant. 

(A,C,E,G,I)  Wild-type.  (B,D,F,H,J)  nflaZFN-mut2;nflbZFN~mut2  double  homozygous  mutant.  (A-D) 
Brightfield  images.  (E-H)  Epi-fluorescent  images.  (I,J)  Confocal  images.  (A,B)  Overall  morphology  of 
the  mutant  looks  normal  at  6  dpf.  (C,D)  Dorsal  view  of  the  tail  region.  Arrows  indicate  missing 
melanocytes  in  the  lateral  strip  of  the  mutant.  (E,F)  nfla;nflb  mutant  larva  has  more  dorsally  migrated 
o/zg2:EGFP+  OPCs  than  control  (bracket).  (G,F)  At  7  dpf,  dbh:G¥V+  SCG  (arrowhead)  is  enlarged  in  the 
mutant  larva.  (I,J)  Transverse  sections,  dorsal  up.  Phospho-ERKl/2  is  highly  detected  in  a  few  neurons 
and  white  matter  of  the  spinal  cord  of  the  nfla;nflb  mutant  larva. 


Page  27 


